Dirac spectra and edge states in honeycomb plasmonic lattices 
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We study theoretically the dispersion of plasmonic honeycomb lattices and find Dirac spectra 
for both dipole and quadrupole modes. Zigzag edge states derived from Dirac points are found 
in ribbons of these honeycomb plasmonic lattices. The zigzag edge states for out-of-plane dipole 
modes are closely analogous to the electronic ones in graphene nanoribbons. The edge states for 
in-plane dipole modes and quadrupole modes, however, have rather unique characters due to the 
vector nature of the plasmonic excitations. The conditions for the existence of plasmonic edge states 
are derived analytically. 
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Metal nanoparticles support remarkable optical reso- 
nances known as localized surface plasmons 0] . Sur- 
face plasmons are collective electronic density waves 
bound at metal surfaces that can locally enhance inci- 
dent optical fields by orders of magnitude. The enhanced 
optical fields in the near field give rise to many extraordi- 
nary opticalphenomena such as surface-enhanced Raman 
scattering [s, Q] and single-molecule fluorescence detec- 
tion The arrangement of metal nanoparticles regu- 
larly can offer more. It has been shown that a chain of 
closely spaced metal nanoparticles can guide electromag- 
netic (EM) energy with lateral mode confinement below 
the diffraction limit P|. Two-dimensional (2D) lattices 
consisting of metal nanoparticles show interesting tun- 
able optical response over the entire visible range Q. 
A three-dimensional arrangement of metal nanoparticles 
can form metallic photonic crystals (PCs) with robust 
photonic band gaps that depend on the local order rather 
than on the symmetry or the global long range order 
Q. Metal nanoparticles thus manifest promising build- 
ing blocks for plasmonic materials, leading to unprece- 
dented applications in nanophotonics, nonlinear optics, 
bio-sensing, and even medical therapy. 

Recently, graphene has received considerable interest 
because of the existence of the Dirac point that offers 
remarkable electronic properties BUS El- The Dirac 
spectra in 2D PCs for EM waves have also been identi- 
fied Ij, |l3[ and the Dirac-point-derived edge states have 
novel transportation properties when time-reversal sym- 
metry breaking is introduced. These edge states can ren- 
der "one-way" waveguiding for EM waves, confirmed by 
numerical simulations 14]. Both graphene and 2D PCs 
support Dirac dispersions. However, the Dirac spectrum 
in graphene is derived directly from the nearest-neighbor 
hopping of bound electron states, while that in 2D PCs 
corresponds to the photonic bands of scattering photons. 

In this paper, we study theoretically the plasmonic 
band structures of metal nanospheres arranged in a hon- 
eycomb lattice. The plasmonic lattice is an open system 
that supports both guided modes and leaky modes [3], 



and as such, it is different from electronic graphene and 
2D PCs. Furthermore, the Dirac dispersions found in 
graphene and 2D PCs are only for scalar waves. The 
2D plasmonic lattice can support both scalar waves and 
vector waves. Our results reveal that for both scalar 
and vector waves, this honeycomb plasmonic lattice pos- 
sesses Dirac points for the infinite system and guided edge 
states when an edge is present. The existence of Dirac 
points and edge states in plasmonic systems may open 
up new avenues in both physics and applications for the 
emerging field of plasmonics. 

The honeycomb plasmonic structure under study is 
shown in Fig. 1(a), in which each black dot in the 
X — y plane represents a metal sphere. Its Brillouin zone 
(BZ) is shown Fig. 1(b). The inter-distance between 
adjacent spheres is uq — a/-\/3 , where a is the lattice 
constant. The primitive and reciprocal lattice vectors 
are, respectively, given by ai_2 — a(±l/2, \/3/2) and 
hi 2 = 27r(±l, l/V3)/a. We use the Drude-type per- 
mittivity e{uj) = 1 — ijJp/i^^ with LUp = 6.18 eV for metal 
spheres. The lattice constant is chosen to be a = 60 
nm, and the sphere radius is = 10 nm. For this ratio 
Vs/a = \/3/6, the dipole and quadruple dispersions are 
separated in energy as demonstrated below. 

The plasmonic dipole dispersions calculated within the 
quasi-static approximation (QSA) limit are shown in Fig. 
1(c), in which the red dots are for out-of-plane Pout 
modes and the black dots correspond to in-plane P,;„ 
modes. The Pi„ and Font modes form their own orthog- 
onal subspaces. We see that there are two Dirac cones at 
the K point, one for the Pm and one for the Pout modes. 
Fig. 1(d) shows the plasmonic dispersions in the QSA 
limit if only quadrupoles (Q) are considered. We found 
a Dirac point at about 3.9 eV, which is the single sphere 
quadruple resonance frequency. At these Dirac points, 
the constant frequency surfaces are represented by two 
cones meeting at the K point. It should be pointed out 
that the in-plane dipole and quadrupole modes are vec- 
tor waves in nature. The existence of Dirac dispersions in 
these cases is a result of the symmetry of the honeycomb 
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FIG. 1: (Color online) (a) Schematic structure of the honey- 
comb plasmonic lattice, ai, a2 are the primitive vectors; a 
and z are the ID primitive vectors for armchair and zigzag 
ribbons respectively, (b) Reciprocal lattice vectors, (c) QSA 
band structures for in-plane dipole modes (black dots) and 
out-of-plane dipole modes (red dots); (d) QSA band struc- 
tures for quadrupoles. MST results (electrodynamic calcula- 
tion with retardation effects taken fully into account) are also 
shown in (c) and (d) by green dots. The dashed lines in pan- 
els (c) and (d) indicate the light lines. The lattice constant is 
60 nm and the sphere radius is 10 nm. 



structure which is robust against their vector nature. 

We note that the l-th multipole plasmon reso- 
nance frequency of a small plasmonic sphere is uji = 

1/2 

[^/(2Z-|-1)] ' LUp. When the spheres are arranged in a 
lattice, each multipole resonance couples to form a set 
of bands, and these manifolds do not mix when Ts/oo is 
small. The dipole and quadrupolc bands shown in Fig. 1 
are indeed separated in energy. In Figs. 1(c) and (d), we 
compare the QSA results with the fully fledged multiple 
scattering theory (MST) results (shown in green color) 
[l^ . The MST gives a solution to Maxwell equations that 
takes into account the full retardation effects. The com- 
parison between the MST (green dots) and QSA results 
is shown in Fig. 1(c) for dipole modes and in Fig. 1(d) 
for quadrupolc modes. The details of the MST method 
are described elsewhere [l^. The MST dispersions are 
shown only for k > lj/c, since the leaky modes inside 
the light cone are not well deflned due to the coupling 
with free photons. As the K point is far away from light 
cone, the Dirac dispersion found by QSA agrees well with 
the MST calculations. We will discuss the Dirac disper- 
sions with QSA since it can give a more straightforward 
interpretation. 

An important consequence of the Dirac spectrum in 
graphene electronics is the existence of peculiar edge 
states 111 with unusual physical properties jl3| in 
graphene ribbons with finite widths. Two types of 
graphene nanoribbons, namely armchair and zigzag rib- 
bons, are usually considered. In Fig. 1(a), the vector 
a (z) indicates the translation vector of the armchair 
(zigzag) ribbons, and the shorter (longer) side of the rect- 




FIG. 2: (Color online) Band structures for TV = 20 zigzag 
plasmonic ribbons (blue) and projected bands for the 2D 
honeycomb structure along the corresponding reciprocal axis 
(grey area), (a) Out-of-plane dipole modes. The flat band 
between 27r/3 and n corresponds to edge states, (b) In- 
plane dipole modes. The flat band with ka ranging from 
to 27r/3 corresponds to Dirac-point-derived edge states, (c) 
Qudrupole modes. Near the single particle resonance 3.9 eV 
(<^i=2), there are two additional bands: one with ka rang- 
ing from to 27r/3, the other traversing the whole BZ. The 
system parameters are the same as in Fig. [T] 



angle inside the 2D graphene BZ is the ID BZ of the arm- 
chair (zigzag) ribbons. The band structure of the arm- 
chair (zigzag) graphene ribbons can be predicted by the 
projection of 2D graphene bands onto the corresponding 
axis of the reciprocal lattice vector. The original Dirac 
point K {K') is expected to appear at fca = (27r/3) for 
armchair (zigzag) ribbons. Compared with the projected 
2D electronic graphene bands, a new feature of zigzag rib- 
bons is the appearance of edge states that extends from 
fca = TT to 27r/3, where fca = 27r/3 is the expected degen- 
eracy point of the 2D projected graphene bands. These 
new additional edge states form a flat band at energy 
E ^ with wave functions localized near the edge sites. 
Within a tight-binding model, the necessary condition 
for _E = edge states is " zero-hopping" , namely, the to- 
tal sum of hopping terms over the nearest neighbor sites 
should vanish. For a semi-inflnite graphene sheet with 
a zigzag edge, the tight-binding model shows that when 
fca = TT, the edge state wave function is completely lo- 
calized. Thereafter the decay length becomes longer and 
smoothly connects to a delocalized state at fca = 27r/3, 
which is the 2D graphene Dirac point. 

We now discuss the Dirac-point-derived edge states 
of plasmonic ribbons. Unless otherwise stated, the re- 
sults for finite-width plasmonic ribbons here are calcu- 
lated for = 20 zigzag and armchair ribbons, where N 
is an integer characterizing the width of the ribbon [ll[ . 
We expect that the edge states of out-of-plane modes 
might be similar to those of electronic edge states de- 
rived from the out-of-plane it states in graphene, while 
the edge states for the Pin and Q modes (if there are 
any) might be very different. We note that electronic 
graphene jll| and 2D PCs studied recently [l^ are all 
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FIG. 3: (Color online) Same as Fig. 2 but for armchair plas- 
monic ribbons. There is no edge state for the armchair ribbon 
near the single particle resonance at 3.52 eV (dipole) or 3.9 
eV (quadrupole) . 



dealing with scalar waves. In their cases, the electron 
and photon transport near the Dirac points can be de- 
scribed by the "two-component" massless Dirac Hamil- 
tonian H = "/{kxCx + kyay), where 7 is a band parameter 
and (7 is the Pauli matrix. In the plasmonic sphere arrays, 
the vector nature of EM waves enters explicitly, and the 
underlying physics can go beyond the "two-component" 
Dirac equation. 

In Fig. 2, the blue dots represent the plasmonic bands 
of zigzag ribbons for out-of-plane dipole, in-plane dipole 
and quadrupole. The bands are superimposed on top of 
the plasmonic bands of the 2D honeycomb lattice pro- 
jected along the zigzag axis, which are shown as the grey 
background. For Font modes in Fig. 2(a), we see that the 
ribbon bands basically overlap the projected 2D bands 
with the exception of an additional degenerate flat band 
starting from about ka — 27r/3 to tt in the gap of the 
projected bands, which corresponds to the edge states in 
the zigzag graphene ribbons. This flat band occurs at 
3.51 eV (the single sphere dipole resonance uji=i — 3.52 
eV), which has the similar physical meaning oi E — 
("zero-hopping") in graphene. The slight differences in 
the edge state frequencies and (jJ;=i come from the hop- 
ping among farther sites in the lattice, which is neglected 
in the tight-binding model for graphene. 

For Pi„ modes, there is no edge state from ka = 27r/3 
to TT, however, a new additional flat band with its fre- 
quency near LUi^i emerges from fca = to about 27r/3. 
This phenomenon is rather different from that in the 
zigzag graphene ribbon, where the edge states are for- 
bidden near the F point. The "two-component" effective 
Dirac equation commonly used in graphene electronics is 
not valid in this case since the vector nature of EM waves 
introduces extra degrees of freedom. More interesting ob- 
servation can be made when we come to qudrupoles. Two 
sets of edge states appear near the single particle reso- 
nance at uji=2 = 3.9 eV. The dispersion of one set of edge 
states is similar to the in-plane dipole edge mode dis- 
cussed above, as it is almost flat and is confined between 



< ka < 27r/3. There is an additional band of edge 
states with a small and negative group velocity that goes 
all the way from fca = to tt. 

The Pout, Pin and Q modes for the armchair plas- 
monic ribbons are shown in Fig. 3 as blue dots, overlaid 
on top of the 2D plasmonic bands (grey color) projected 
along the armchair axis. There is no additional flat band 
near the single sphere resonance a;/=i for dipoles or w/=2 
for quadrupoles. We can thus conclude that the Dirac- 
point-derived edge states observed in the plasmonic sys- 
tem come from the boundary condition (zigzag edge) 
which is analogous to that in electronic graphene. 

In the following, we give an analytic derivation of 
the condition for the existence of plasmonic zigzag edge 
states near multipole resonance frequencies by truncating 
the Green function to the nearest neighbor. Let us con- 
sider a general tensor field F of rank I (which can be the 
electric field in the case of a dipole or the field gradient in 
the case of a quadrupole) generated by a localized source 

V with rank I (dipole, quadrupole . . . ). In free space, the 
field and the source are related by the 2^-th rank Green 
function G(r — r') which satisfies F(r) = G(r — r')V(r') 
and the linear response gives V(r) = aF(r), where a is 
the multipole polarizability. We place these point sources 

V on a semi-infinite honeycomb lattice with a zigzag edge 
as shown in the inset of Fig. 4. The indices (to, n) in 
the inset are used to label the honeycomb sites, where 
n and to label the position along (or y-axis) and per- 
pendicular to (x-axis) the zigzag chain respectively (see 
Fig. 4). rm,7i is the coordinate of the (m, n) site, and the 
site (1, n -I- 1/2), located between (1, n) and (1, n + 1), is 
chosen to be the origin for convenience. We search for 
eigenstates with a Bloch wavevector k along the zigzag 
edge and we demand that the frequency should be equal 
to the single particle multipole resonance frequency. If 
we truncate the Green function to the nearest neighbor, 
the condition for the frequency to be equal to the single 
particle multipole resonance frequency is " zero- hopping" , 
namely a zero sum of the fields from the three nearest 
neighbors. This gives 

G(-ri,„)Vi,„+G(-ri,„+i)Vi,„+i+G(-r2.„+i)V2,„+i = 0, 

(1) 

where \'m,n ~ V,„e™'^" if the periodic Bloch condition 
is imposed. From Eq. (1) a matrix T relating Vi and 
V2 can be defined as 

T - -G-i(-r2.„+i)[G(-ri.„)e-''='^/2+G(-ri,„+i)e''='^/2]. 

(2) 

This procedure can be done iteratively to obtain = 
rpm-i rpj^g eigenvalues of the matrix T give the prop- 
erties of the states. If |A| < 1, the states are localized 
edge states, while |A| = 1 corresponds to a Bloch state 
extending to the interior of the ribbon. 

For out-of-plane modes, the Green function G^°"* (r — 
r') is just an isotropic scalar function and Eq. (1) is 
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FIG. 4: (Color online) Absolute values of eigenvalues A of the 
matrix T as a function of the Bloch wavevector k. V is ei- 
ther P or Q. |A| < 1 gives the condition for edge states, (a) 
Dipole modes. The green line is for out-of-plane modes and 
the red and black lines are for in-plane modes, (b) Quadrupole 
modes. Qi (blue) is always an edge state, while Pi and Q2 
give edge state solutions when ka < 2n/3. The inset is a 
sketch of a semi-infinite honeycomb structure with a zigzag 
edge. The indices (m, n) and vector wave functions are de- 
fined as shown. 



reduced to p.ie'"'"' + p^^(,Kn+i)ka 



p^2e'^""' 2/-- = 0, 
which is exactly the same condition found for electronic 
graphene ll|. This is why the edge state behavior of 
the out-of-plane mode is so similar to that of electronic 
graphene ribbons. The ratio PzijPzX is simply given 
by — 2cos/ca and is shown by green lines in Fig. 4(a). 
For the in-plane dipole mode, the QSA dyadic dipolar 
Green function can be written as GP'"(r — r')Pi„(r') = 
[3(n • Pj„)n — Pi„] / |r — r'l , where n is the unit vec- 
tor along the r — r' direction, and T is a 2 x 2 matrix 
that can be found analytically. In Fig. 4(a), the ab- 
solute values of the two eigenvalues are plotted by red 
and black lines. The black curve is always larger than 1, 
while the red curve is smaller than 1 when ka < 2n/3, 
denoted as Pi in the figure, which explains why in-plane 
dipolar edge states are confined by < ka < 27r/3. 
For quadrupoles, there are five independent variables 
{Qyy,Qzz,Qxy,Qxz,Qyz} and the Corresponding QSA 
Green function can be derived analytically after tedious 
algebra. The details will be published elsewhere. T- 
matrix becomes a 5 x 5 matrix, and the absolute values 
of the five eigenvalues are plotted in Fig. 4(b). Only two 
sets of eigenvalues give absolute values less than 1 and 
thus form edge states. The edge state solution labeled 
as Qi (blue) is smaller than unity in the whole ID BZ. 
This corroborates the existence of an edge state band 
that goes from fca = to tt in Fig. 2(c) in the TV = 20 
zigzag plasmonic ribbon. Another state denoted by Q2 
(red) has a critical value for the edge state at ka = 27r/3 
similar to the in-plane dipole mode Pi. 

Another interesting observation that can be deduced 
from Fig. 4 is that the eigenvalues for the T-matrix go 
to zero at ka — 0.267r and 0.547r for Pi and Qi respec- 
tively, which implies that the corresponding edge states 
are localized at the edge completely. In contrast, com- 



pletely localized states in electronic graphene are only 
found at the BZ boundary ka = n. The completely lo- 
calized state away from the BZ boundary comes from the 
vector nature of EM waves. 

In summary, we have found Dirac spectra in 
"graphene-like" plasmonic structures for both dipole and 
quadrupole coupling. Edge states for plasmonic lattices 
with finite widths are also identified and analyzed numer- 
ically and their existence conditions are derived analyti- 
cally. The edge states for the out-of-plane dipole modes 
in zigzag plasmonic ribbons resemble those of the elec- 
tronic grapheme zigzag edge states, but the edge states 
for in-plane dipole and quadrupole modes have rather 
unique characters due to their vector nature. These re- 
sults can lead to the study of Dirac point behavior in yet 
another class of interesting material. 
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